Quantifying defects in graphene via Raman spectroscopy at different excitation 

energies 
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We present a Raman study of Ar+-bombarded graphene samples with increasing ion doses. This 
allows us to have a controlled, increasing, amount of defects. We find that the ratio between the 
D and G peak intensities for a given defect density strongly depends on the laser excitation energy. 
We quantify this effect and present a simple equation for the determination of the point defect 
density in graphene via Raman spectroscopy for any visible excitation energy. We note that, for 
all excitations, the D to G intensity ratio reaches a maximum for an inter-defect distance ~ 3 nm. 
Thus, a given ratio could correspond to two different defect densities, above or below the maximum. 
The analysis of the G peak width and its dispersion with excitation energy solves this ambiguity. 



I. INTRODUCTION 



Quantifying defects in graphene related systems, which 
include a large family of sp^ carbon structures, is cru- 
cial both to gain insight in their fundamental proper- 
ties, and for applications. In graphene, this is a key 
step towards the understanding of the limits to its ul- 
timate mobilityi"—. Large efforts have been devoted 
to quantify defects and disorder using Raman spec- 
troscopy for nanographites^"— , amorphous carbonsi^r— , 
carbon nanotubes^"*'^^, and graphen o^^'^^" — . The first 
attempt was the pioneering work of Tuinstra and Koenig 
(TK)^. They reported the Raman spectrum of graphite 
and nano-crystalline graphite, and assigned the mode at 
~1580cm~^ to the high frequency i?2g Raman allowed 
optical phonon, now known as G peak^. In defected and 
nanocrystalline samples they measured a second peak at 
^1350cm^^, now known as D peak-. They assigned it 
to an Aig breathing mode at the Brillouin Zone (BZ) 
boundary K, activated by the relaxation of the Raman 
fundamental selection rule q « 0, where q is the phonon 
wavevector— . They noted that the ratio of the D to G 
intensities varied inversely with the crystallite size, ia- 
Refjii noted the failure of the TK relation for high de- 
fect densities, and proposed a more complete amorphiza- 
tion trajectory valid to date. Refsi^'^ i^''''^^ reported a 
significant excitation energy dependence of the intensity 
ratio. Refsi^iiS measured this excitation laser energy de- 
pendency in the Raman spectra of nanographites, and 
the ratio between the D and G bands was shown to de- 
pend on the fourth power of the excitation laser energy 

There is, however, a fundamental geometric difference 
between defects related to the size of a nano-crystallite 
and point defects in the sp^ carbon lattices, resulting 
in a different intensity ratio dependence on the amount 
of disorder. Basically, the amount of disorder in a 



nano-crystallite is given by the amount of border (one- 
dimensional defects) with respect to the total crystallite 
area, and this is a measure of the nano-crystallite size La- 
in graphene with zero-dimensional point-like defects, the 
distance between defects, Ld, is a measure of the amount 
of disorder, and recent experiments show that different 
approaches must be used to quantify Ld and La by Ra- 
man spectroscopy^!. The effect of changing Ld on peak 
width, frequency, intensity, and integrated area for many 
Raman peaks in single layer graphene was studied in 
Refj^, and extended to N-layer graphene in Refi^, all 
using a single laser line = 2.41 eV. 

Here, to fully accomplish the protocol for quantifying 
point-like defects in graphene using Raman spectroscopy 
(or equivalently, Ld), we use different excitation laser 
lines in ion-bombarded samples and measure the D to G 
peak intensity ratio. This ratio is denoted in literature as 
Iu/Ig or I(D)/I(G), while the ratio of their areas, i.e. fre- 
quency integrated intensity, as Au/Aq or A(D)/A(G). In 
principle, for small disorder or perturbations, one should 
always consider the area ratio, since the area under each 
peak represents the probability of the whole process, con- 
sidering uncertainty^ii^^. However, for large disorder, it 
is far more informative to decouple the information on 
peak intensity and full width at half maximum. The lat- 
ter, denoted in literature as FWHM or F, is a measure 
of structural disorder-i2i2ii^, while the intensity repre- 
sents the phonon modes/molecular vibrations involved in 
the most resonant Raman processe s^^i^^'^^ . For this rea- 
son, in this paper we will consider the decoupled Id / Ig 
and peak widths trends. We find that, for a given Ld, 
Id/Iq increases as the excitation laser energy increases. 
We present a set of empirical formulas that can be used 
to quantify the amount of point-like defects in graphene 
samples with Ld > 10 nm using any excitation laser en- 
ergy/wavelength in the visible range. The analysis of the 
D and G peak widths and their dispersions with excita- 
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FIG. 1. Raman spectra of five ion bombarded SLG measured 
at El, =2.41eV (Al = 514.5 nm). The Lu values are given 
according to Ref.— , and the main peaks are labeled. The 
notation within parenthesis [e.g. 2D(G')] indicate two com- 
monly used notations for the same peak (2D and G' )^°'^° . 



tion energy unambiguously discriminate between the two 
main stages of disordering incurred by such samples. 



II. RESULTS AND DISCUSSION 

We produce single layer graphene (SLG) samples with 
increasing defect density by mechanical exfoliation fol- 
lowed by Ar^-bombardment, as for the procedure out- 
lined in Refi^. The ion-bombardment experiments are 
carried out in an OMICRON VT-STM ultra-high vac- 
uum system (base pressure 5xlO~^^mbar) equipped 
with an ISE 5 Ion Source. Raman spectra are mea- 
sured at room temperature with a Renishaw micro- 
spectrometer. The spot size is ^ l/im for a 100 x objec- 
tive, and the power is kept at ^ 1.0 mW to avoid heating. 
The excitation energies, i^L, (wavelengths, Al) are: Ti- 
Sapph 1.58eV (785 nm), He-Ne 1.96eV (632.8nm), Ar+ 
2.41 eV (514.5 nm). 

Figure [1] plots the Raman spectra of five SLG ex- 
posed to different ion bombardment doses in the range 
10" Ar+/cm2 (one defect per 4 x lO" C atoms) to 10^^ 
Ar+/cm^ (one defect for every four C atoms). The bom- 
bardment procedure described in Ref.'^^ is accurately 
reproducible. By tuning the bombardment exposure 
we generated samples with Ld — 24, 14, 13, 7, 5, and 2 nm. 
All spectra in Fig. [T]are taken at i?L =:2.41eV (Al = 
514.5 nm). 

The Raman spectra in Figure [T] consist of a set of dis- 



tinct peaks. The G and D appear around 1580 cm~^ and 
1350 cm~^, respectively. The G peak corresponds to the 
i?2g phonon at the Brillouin zone center. The D peak is 
due to the breathing modes of six-atom rings and requires 
a defect for its activation^'iiii^'^. It comes from trans- 
verse optical (TO) phonons around the K or K' points in 
the 1*** Brillouin zone^ i^'''^^ , involves an intervalley double 
resonance process^i^, and is strongly dispersive'^* with 
excitation energy due to a Kohn Anomaly at K"^^ . Dou- 
ble resonance can also happen as intravalley process, i. e. 
connecting two points belonging to the same cone around 
K or K'"^^. This gives the so-called D' peak, which is cen- 
tered at ~ 1620 cm~^ in defected samples measured at 
514.5nroi^. The 2D peak (also called G' in the litera- 
ture) is the second order of the D peal i^^'^° . This is a 
single peak in single layer graphene, whereas it splits in 
four in bilayer graphene, reflecting the evolution of the 
electron band structure^SdS. The 2D' peak (also called 
G" in analogy to G') is the second order of D'. Since 
2D(G') and 2D'(G") originate from a process where mo- 
mentum conservation is satisfied by two phonons with 
opposite wavevectors, no defects are required for their 
activation, and are thus always present. On the other 
hand, the D -|-D' band (^ 2940 cm~^) is the combination 
of phonons with different momenta, around K and F, 
thus requires a defect for its activation. 

Refill proposed a three stage classification of disorder 
in carbon materials, to simply assess the Raman spec- 
tra of carbons along an amorphization trajectory lead- 
ing from graphite to tetrahedral amorphous carbon: 1) 
graphite to nanocrystalline graphite; 2) nanocrystalline 
graphite to low sp'^ amorphous carbon; 3) low sp'^ amor- 
phous carbon to high sp^ (tetrahedral) amorphous car- 
bon. In the study of graphene, stages 1 and 2 are the 
most relevant and are summarized here. 

In stage 1, the Raman spectrum evolves as fol- 
lowsiii^i^^: a) D appears and Id/Ig increases; b) D' 
appears; c) all peaks broaden. In the case of graphite 
the D and 2D lose their doublet structureili^; e) D D' 
appears; f) at the end of stage 1, G and D' are so wide 
that they start to overlap. If a single lorentzian is used 
to fit G -f D', this results in an upshifted wide G band at 
-1600cm-i. 

In stage 2, the Raman spectrum evolves as foUowsii: 

a) the G peak position, denoted in literature as Pos(G) 
or ujQ , decreases from 1 600 cm~ ^ towards ~ 1 5 1 cm~ ^ ; 

b) the TK relation fails and Iu/Ig decreases towards 0; 

c) u)Q becomes dispersive with the excitation laser en- 
ergy, the dispersion increasing with disorder; d) there 
are no more well defined second-order peaks, but a small 
modulated bump from ^ 2300 cm^^ to ^ 3200cm~^— i^^. 

In disordered carbons ujq increases as the excitation 
wavelength decreases, from IR to UV— . The disper- 
sion rate, Disp(G) = Acjc/Ai^L, increases with disor- 
der. The G dispersion separates the materials into two 
types. In those with only sp^ rings, Disp(G) saturates at 
^ 1600 cm"-'^, the G position at the end of stage 1. In con- 
trast, for those containing sp^ chains (such as in amor- 
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phous and diamond-like carbons), G continues to rise 
past 1600 cm~^ and can reach ~1690cm~^ for 229 nm 
excitationiiiiS,. On the other hand, D always disperses 
with excitation energ y^^i^^ . Tq always increases with 
disorderi2i2^i22i2^. Thus, combining Id/Iq and Tq al- 
lows to discriminate between stages 1 or 2, since samples 
in stage 1 and 2 could have the same I^/Ig, but not the 
same Tq, being this much bigger in stage 3^i21i^. 

We note that Figure [1] shows the loss of sharp second 
order features in the Raman spectrum obtained from the 
Ld = 2 nm SLG. This is an evidence that the range of 
defect densities in our study covers stage 1 (samples with 
Ld = 24, 14, 13, 7, 5 nm) and the onset of stage 2 (sample 
with Ld = 2 nm) . 

Figures [2^-c report the first-order Raman spectra of 
our ion-bombarded SLGs measured at Ei^ =1.58eV 
(Al =785nm), 1.96eV (632.8nm),2.41 eV (514.5nm), 
respectively. Figure shows the Raman spectra of 
the ion-bombarded SLG with — 7 nm obtained using 
the three different laser energies. We note that Iu/Ig 
considerably changes with the excitation energy. This 
is a well-know effect in the Raman scattering of sp^ 
carbons^iiaiiiii^ii^ii^. RefJ^ noted that the integrated 
areas of different peaks depend differently on excitation 
energy Ei^: while Ad, ^d', and A2Y) shown no i?L- 
dependence, Aq was found to be proportional to E^. The 
independence of A2D on i?L agrees with the theoretical 
predictionii if one assumes that the electronic scatter- 
ing rate is proportional to the energy. However, a fully 
quantitative theory is not trivial since, in general, Ad 
depends not only on the concentration of defects, but 
on their type as well (e.g., only defects able to scatter 
electrons between the two valleys can contribute)^i2^i^. 
Different defects can also produce different frequency and 
polarization dependence of Ajj^^^^i^^^. 

Figure |3] plots Id/Ig for all SLGs and laser ener- 
gies. For all E'l, Id/Ig increases as Xd decreases (stage 
1), reaches a maximum at Xd ~3nm, and decreases to- 
wards zero for iD<3nm (stage 2). It is important 
to understand what the maximum of Id /Iq vs. Ld 
means. Id will keep increasing until the contribution 
from each defect sums independently^li^i. In this regime 
(stage 1) Id is proportional to the total number of de- 
fects probed by the laser spot. For an average defect 
distance Ld and laser spot size Ll, there are on aver- 
age (Ll/^d)^ defects in the area probed by the laser, 
thus /d oc (Ll/^d)^- On the other hand, /q is propor- 
tional to the total area probed by the laseroc Lf^, giving 
/d//goc 1/L p^^'^^ . However, if two defects are closer 
than the average distance an e-h pair travels before scat- 
tering with a phonon, then their contributions will not 
sum independently anymore^L^^i^i^^. This distance can 
be estimated as vf/lud ~3nm^, where vf ~ lO^m/s is 
the Fermi velocity around the K and K' points, in ex- 
cellent agreement with the predictions of Refs.— and the 
data of Refsi ^^'^^1'^^ . For an increasing number of defects 
(stage 2), where Ld < 3nm, sp^ domains become smaller 
and the rings fewer and more distorted, until they open 



up. As the G peak is just related to the relative motion 
of sp^ carbons, we can assume Iq roughly constant as a 
function of disorder. Thus, with the loss of sp^ rings. Id 
will decrease with respect to Iq and the Id //g oc 1 /Lp 
relation will no longer hold. In this regime. Id /Ig oc M 
(M being the number of ordered rings), and the develop- 
ment of a D peak indicates ordering, exactly the opposite 
to stage 1™. This leads to a new relation: /d//goc Ld— . 

The solid lines in Fig. |3] are fitting curves following the 
relation proposed in Ref<21: 



Id 
Ig 



= Ca 



4) 



2 ,t2 



(1) 

The parameters ta and rs are length scales which de- 
termine the region where the D band scattering takes 
place, rs determines the radius of the structurally disor- 
dered area caused by the impact of an ion. ta is defined 
as the radius of the area surrounding the point defect in 
which the D band scattering takes place, although the sp^ 
hexagonal structure is preserved^^. In short, the differ- 
ence TA — rs defines the Raman relaxation length of the 
D band scattering, and is associated with the coherence 
length of electrons which undergo inelastic scattering by 
optical phonons^^''^'^. The fit in Figure[2]is done consider- 
ing rs — l nm (as determined in Refi^ and expected to be 
a structural parameter, i.e. Ei, independent). Further- 
more, within experimental accuracy, all data can be fit 
with the same rA = 3.1nm, in excellent agreement with 
the values obtained in Refsi22i2^iH. Any uncertainty in 
TA does not affect the results in the low defect density 
regime (Ld > 10 nm) discussed later. 

Ref.^^ suggested that Id/Ig depends on both an ac- 
tivated (A) area, pounded by the parameter Ca, and a 
structurally defective area (S), pounded by a parameter 
Cs- Here we selected Cs = in eq. H]) for two reasons: 
(i) Cs should be defect-structure dependent, and in the 
ideal case where the defect is the break-down of the G-G 
bonds, Cs should be null; (ii) here we do not focus on the 
large defect density regime, Ld < rs . The parameter Ca 
in eq. ^ corresponds to the maximum possible Id/Ig, 
which would be observed in the ideal situation where the 
D band would be activated in the entire sample with no 
break down of any hexagonal carbon ringSI. 

Ca has been addressed in Refi^ as related to the ra- 
tio between the scattering efficiency of optical graphene 
phonons evaluated between F and K. As we show here, 
the large Id /Ig dependence on L^l comes from the change 
on Ca, which suggests this parameter might also depend 
on interference effects, when summing the different elec- 
tron/hole scattering processes that are possible when ac- 
counting for the Raman cross section^"—. Note that Ca 
decreases as the laser energy increases. The solid line in 
the inset to Fig.[2]is the fit of the experimental data (dark 
squares) by using an empirical relation between the max- 
imum value of Id/Ig and L^l, of the form Ca = AE^^ . 
The fit yields A = (160 ± 48) eV, by setting 5 = 4 in 
agreement with Refsi^ii^. 



4 






'-□=2nm 




L|^=5nm /, 




A 


/-j^=14nm 




Lp=24nm a 






1200 1300 1400 1500 

Raman shift {cm"^) 



1600 1200 1300 1400 1500 1600 1700 1200 1300 1400 1500 1600 1700 1200 1300 1400 1500 1600 1700 

Raman shift (cm"') Raman shift (cm"') Raman shift (cm"') 



FIG. 2. (a-c) Raman spectra of five distinct ion-bombarded graphene samples using the excitation laser energies (wavelengths) 
El =f.58eV (Al =785nm), =i.96eV (Al = 632.8 nm), and El =2.4ieV (Al =5i4.5nm), respectively, (d) Raman 
spectra of an ion-bombarded sample with Lu = 7nm obtained using these three excitation laser energies. 
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FIG. 3. Id/Ig for all samples and laser energies consid- 
ered here. Solid lines are fits according to equation [T] with 
rs = 1 nm, Cs=0, and rA = 3.Inm. The inset plots C'a 
as a function of E-^. The solid curve is given by Ca = 



(160 ± 48) El 



We no"w focus on the lovif-defect density regime (Xd > 
10 nm), since this is the case of most interest in order to 
understand ho"w Raman active defects limit the ultimate 
mobility of graphene samplesi'"— . In this regime, "where 
Ld > 2rA, the total area contributing to the D band 
scattering is proportional to the number of point defects, 
giving rise to Ib/Ig oc l/i^, as discussed above. For 



large values of Ld, cq. ([T]) can be approximated to 



r2 



(2) 



By taking ta = 3.1 nm, rs = 1 nm, and also the relation 
Ca = (160 ± 48)£'l obtained from the fit of the exper- 
imental data sho"wn in Figure [5J eq. ^ can be re"written 
as 



(nm^ 



(4.3 ± 1.3) X 10^ ^Jd 



(3) 



In terms of excitation laser "wavelength Al (in nanome- 
ters), we have 



Li, (nm") = (1.8 ± 0.5) x lO^'^A^ ( ^ 



(4) 



Equations ([3]) and ^ are valid for Raman data obtained 
from graphene samples "with point defects separated by 
Ld > 10 nm using excitation lines in the visible range. 
In terms of defect density n-o{cm~^) — IQ^'^ / {tiL^), eqs. 
dSl) and (gl) become 



nD(cm 



(7.3 ±2.2) X m'El — 



and 



7iD(cm ^) = 



(1.8 ±0.5) X 10^ 



(5) 



(6) 



Figure m plots E^{Iy) / Iq) as a function of for the 
data sho"wn in Figure [5] The data -with Ld > 10 nm 
obtained "with different laser energies collapse in the same 
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FIG. 4. Ef^{Ir> / Ig) as a function of Ld for the data shown in 
Figure [21 The dashed blue hne is the plot obtained from the 
substitution of the relation Ca ~ (160)/-/?^^ equation [1] 
The solid dark line is the plot of the product E^{Io/Ig) as 
a function of Ld according to equation [S] The shadow area 
accounts for the upper and lower limits given by the ±30% 
experimental error. 



curve. The dashed blue line is the plot obtained from the 
substitution of the relation Ca = (160)/i?L in eq.[TJ The 
solid dark line is the plot Ef^{lY)/lQ) versus Ld according 
to eqs. dS]) and (|4]). The shadow area accounts for the 
upper and lower limits given by the ±30% experimental 
error. The plot in Fig. 2] validates these relations for 
samples with Ld > 10 nm. 

Figure [5^ plots Fd and F2D as a function of Ld. 
Within the experimental error, a dependence of Fd 
or F2D on the excitation energy during stage 1 can 
not be observed. D and 2D always disperse with 
excitation energy, with Awd/AE'l ~ 52 cm~^/eV, and 

Aw2d/al;l = 2Awd/al;l. 

Figures 5b, c plot the G peak dispersion Disp(G) = 
Aa;G/A£;L and Fq =FWHM(G) as a function of Ld, re- 
spectively. As shown in Figure [SJd, Awq/AL^l remains 
zero until the onset of stage two, when it becomes slightly 
dispersive (AclIq/AL^l '^Scm^^/eV). Fq (Figure [S};) re- 
mains roughly constant at ^ 14cm~^, a typical value for 
as-prepared exfoliated graphen o^^i^°i^°'^^ , until the on- 
set of stage 2 (corresponding to the maximum Ld/Lq) 
as suggested in Refi^, and shown in Refi^ for a single 
laser line L^l = 2.41 eV. Combining Iu/Ig and Fq allows 
to discriminate between stages 1 or 2, since samples in 
stage 1 and 2 could have the same I-q/Iq, but not the 
same Fq, which is much larger in stage 2^^i^^ . 
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FIG. 5. (a) Plot of Fd and versus Ld- (b) G peak disper- 
sion [Disp(G) = Acjc/AiJL] as a function of Ld. Ai^g/ASl 
remains zero until the onset of stage 2. (c) FWIIM(G) = Fg 
as a function of Ld. As suggested in Refs .'^'^'^^ , Fg remains 
roughly constant until the onset of the second stage of amor- 
phization, corresponding to the maximum Id/Ig- 



III. CONCLUSIONS comes from the parameter Ca- 



In summary, wc discussed the use of Raman spec- 
troscopy for quantifying the amount of point-hke defects 
in graphene. We used different excitation laser hnes in 
ion-bombarded samples in order to measure their respec- 
tive Id/Ig- We find that Iu/Ig, for a specific Ld, de- 
pends on the laser energy. We presented a set of em- 
pirical relations that can be used to quantify point de- 
fects in graphene samples with Ld > 10 nm via Raman 
spectroscopy using any laser line in the visible range. 
We show that the Raman coherence length ta is Ej,- 
independent, while the strong Ej, dependence for I-d/Iq 
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